Partial differential equation models for invasive species spread in the presence of spatial heterogeneity.

Models of invasive species spread often assume that landscapes are spatially homogeneous; thus simplifying analysis but potentially reducing accuracy. We extend a recently developed partial differential equation model for invasive conifer spread to account for spatial heterogeneity in parameter values and introduce a method to obtain key outputs (e.g. spread rates) from computational simulations. Simulations produce patterns of spatial spread which appear qualitatively similar to observed patterns in grassland ecosystems invaded by exotic conifers, validating our spatially explicit strategy. We find that incorporating spatial variation in different parameters does not significantly affect the evolution of invasions (which are characterised by a long quiescent period followed by rapid evolution towards to a constant rate of invasion) but that distributional assumptions can have a significant impact on the spread rate of invasions. Our work demonstrates that spatial variation in site-suitability or other parameters can have a significant impact on invasions and must be considered when designing models of invasive species spread.

Drought- and heat-induced mortality of conifer trees is explained by leaf and growth legacies.

An increased frequency and severity of droughts and heat waves have resulted in increased tree mortality and forest dieback across the world, but underlying mechanisms are poorly understood. We used a common garden experiment with 20 conifer tree species to quantify mortality after three consecutive hot, dry summers and tested whether mortality could be explained by putative underlying mechanisms, such as stem hydraulics and legacies affected by leaf life span and stem growth responses to previous droughts. Mortality varied from 0 to 79% across species and was not affected by hydraulic traits. Mortality increased with species' leaf life span probably because leaf damage caused crown dieback and contributed to carbon depletion and bark beetle damage. Mortality also increased with lower growth resilience, which may exacerbate the contribution of carbon depletion and bark beetle sensitivity to tree mortality. Our study highlights how ecological legacies at different time scales can explain tree mortality in response to hot, dry periods and climate change.

Hybridization mediated range expansion and climate change resilience in two keystone tree species of boreal forests.

Current global climate change is expected to affect biodiversity negatively at all scales leading to mass biodiversity loss. Many studies have shown that the distribution of allele frequencies across a species' range is often influenced by specific genetic loci associated with local environmental variables. This association reflects local adaptation and allele changes at those loci could thereby contribute to the evolutionary response to climate change. However, predicting how species will adapt to climate change from this type of data alone remains challenging. In the present study, we combined exome capture sequences and environmental niche reconstruction, to test multiple methods for assessing local adaptation and climate resilience in two widely distributed conifers, Norway spruce and Siberian spruce. Both species are keystone species of the boreal forest and share a vast hybrid zone. We show that local adaptation in conifers can be detected through allele frequency variation, population-level ecological preferences, and historical niche movement. Moreover, we integrated genetic and ecological information into genetic offset predictive models to show that hybridization plays a central role in expanding the niche breadth of the two conifer species and may help both species to cope better with future changing climates. This joint genetic and ecological analysis also identified spruce populations that are at risk under current climate change.

The rhizosphere microbiome reduces the uptake of arsenic and tungsten by Blechnum orientale by increasing nutrient cycling in historical tungsten mining area soils.

The growth of pioneer plants in metal mining area soil is closely related to their minimal uptake of toxic elements. Pioneer plants can inhibit the uptake of toxic elements by increasing nutrient uptake. However, few studies have focused on the mechanisms by which the rhizosphere microbiome affect nutrient cycling and their impact on the uptake of toxic elements by pioneer plants. In this study, we selected Blechnum orientale to investigate the potential roles of the rhizosphere microbiome in nutrient cycling and plant growth in a historical tungsten (W) mining area. Our results showed that while the arsenic (As) and W contents in the soil were relatively high, the enrichment levels of As and W in the B. orientale were relatively low. Furthermore, we found that the As and W contents in plants were significantly negatively correlated with soil nutrients (S, P and Mo), suggesting that elevated levels of these soil nutrients could inhibit As and W uptake by B. orientale. Importantly, we found that these nutrients were also identified as the most important factors shaping rhizosphere microbial attributes, including microbial diversity, ecological clusters, and keystone OTUs. Moreover, the genera, keystone taxa and microbial functional genes enriched in the rhizosphere soils from mining areas played a key role in nutrient (S, P and Mo) bioavailability, which could further increase the nutrient uptake by B. orientale. Taken together, our results suggest that rhizosphere microorganisms can improve pioneer plant growth by inhibiting toxic element accumulation via the increase in nutrient cycling in former W mining areas.

Responses of vascular plant fine roots and associated microbial communities to whole-ecosystem warming and elevated CO2 in northern peatlands.

Warming and elevated CO2 (eCO2) are expected to facilitate vascular plant encroachment in peatlands. The rhizosphere, where microbial activity is fueled by root turnover and exudates, plays a crucial role in biogeochemical cycling, and will likely at least partially dictate the response of the belowground carbon cycle to climate changes. We leveraged the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experiment, to explore the effects of a whole-ecosystem warming gradient (+0°C to 9°C) and eCO2 on vascular plant fine roots and their associated microbes. We combined trait-based approaches with the profiling of fungal and prokaryote communities in plant roots and rhizospheres, through amplicon sequencing. Warming promoted self-reliance for resource uptake in trees and shrubs, while saprophytic fungi and putative chemoorganoheterotrophic bacteria utilizing plant-derived carbon substrates were favored in the root zone. Conversely, eCO2 promoted associations between trees and ectomycorrhizal fungi. Trees mostly associated with short-distance exploration-type fungi that preferentially use labile soil N. Additionally, eCO2 decreased the relative abundance of saprotrophs in tree roots. Our results indicate that plant fine-root trait variation is a crucial mechanism through which vascular plants in peatlands respond to climate change via their influence on microbial communities that regulate biogeochemical cycles.

To what extent can we predict variation of bryophyte and tracheophyte community composition at fine spatial scale along an elevation gradient?

Mounting evidence points to the need for high-resolution climatic data in biodiversity analyses under global change. As we move to finer resolution, other factors than climate, including other abiotic variables and biotic interactions play, however, an increasing role, raising the question of our ability to predict community composition at fine scales. Focusing on two lineages of land plants, bryophytes and tracheophytes, we determine the relative contribution of climatic, non-climatic environmental drivers, spatial effects, community architecture and composition of one lineage to predict community composition of the other lineage, and how our ability to predict community composition varies along an elevation gradient. The relationship between community composition of one lineage and 68 environmental variables at 2-25 m spatial resolution, architecture and composition of the other lineage, and spatial factors, was investigated by hierarchical and variance partitioning across 413 2x2m plots in the Swiss Alps. Climatic data, although significant, contributed less to the model than any other variable considered. Community composition of one lineage, reflecting both direct interactions and unmeasured (hidden) abiotic factors, was the best predictor of community composition of the other lineage. Total explained variance substantially varied with elevation, underlining the fact that the strength of the species composition-environment relationship varies depending on environmental conditions. Total variance explained increased towards high elevation up to 50 %, with an increasing importance of spatial effects and vegetation architecture, pointing to increasing positive interactions and aggregated species distribution patterns in alpine environments. In tracheophytes, an increase of the contribution of non-climatic environmental factors was also observed at high elevation, in line with the hypothesis of a stronger environmental control under harsher conditions. Further improvements of our ability to predict changes in plant community composition may involve the implementation of historical variables and higher-resolution climatic data to better describe the microhabitat conditions actually experienced by organisms.

Understanding boron toxicity in aquatic plants (Salvinia natans and Lemna minor) in the presence and absence of EDTA.

Even though boron is a widely used element in various industries and a contributor to water pollution worldwide, few studies have examined the toxicity of boron in aquatic plants. EDTA is used to maintain aquatic plants cultures, however it is possible to modify the toxicity of metals. The objective of this study is to assess the toxicity of boron in aquatic plants and explore the impact of EDTA presence on the resulting toxic responses. Floating watermoss Salvinia natans and duckweed Lemna minor were exposed to concentrations ranging from 5 to 100 mg/L for 7 days and 1 to 60 mg/L for 3 days, respectively. Growth and photosynthetic activity parameters were investigated in the presence and absence of EDTA. Growth inhibitions in both aquatic plants were observed in a concentration-dependent manner, irrespective of the presence or absence of EDTA. For instance, based on the specific growth rate (leaves coverage), EC10 values for S. natans were calculated as 12.7 (9.9-15.3) mg/L and 8.0 (5.8-10.3) mg/L with and without EDTA, respectively. In the case of L. minor, EC10 values were calculated as 1.3 (0.8-1.89) mg/L and 2.0 (0.4-4.3) mg/L with EDTA without EDTA, respectively. Significant effects were also observed on the photosynthetic capacity, however there was no change in the increase of boron concentration. Generally, negligible effects of EDTA to the toxicity of boron were observed in the present study. By comparing toxicity results based on the presence and absence of EDTA, which is an essential element in the test medium, the results of this study are expected to be utilized for the ecological risk assessment of boron in aquatic ecosystems.

Exploring the effect of different typical plant community on human stress reduction: a field experiment.

Research has demonstrated the positive effect of natural environment on human restoration and well-being. Time spent in nature can often alleviate both physiological and psychological stress. However, few studies have discussed the environmental health effects of the nature's components and characteristics. Sixty volunteers were recruited and one manufactured environment and five different natural environments were randomly assigned to them, including coniferous forests (pure coniferous forest-PC and mixed coniferous forest-MC), broad-leaved forests (pure broad-leaved forest-PB and mixed broad-leaved forest-MB), and mixed forest (mixed coniferous and broad-leaved forest-MCB). Each volunteer sat in a built or natural environment and looked around the environment for 15 min. Physiological (HR, HRV, BP, pulse rate and salivary cortisol) and psychological indicators (POMS and STAI) were used to evaluate the changes in their stress level. Results indicated a strong difference in HR, HRV, POMS and STAI between the built and natural environment, which showed that natural environment can lower the stress level. MC had the best effect on relieving physiological stress, whereas MCB is most successful in improving emotional state and reducing anxiety. Broad-leaved forest and mixed forest significantly affected the DBP and vigor level of the subjects, respectively. While coniferous forest did significantly increase the concentration of salivary cortisol in subjects. The study confirmed that compared to the built environment, the natural environment can relieve the human body's physical and psychological stress and negative emotions, while significantly increasing vitality. And different plant communities also have different effects on the physiological and psychological indicators of the subjects. These results will provide scientific basis for the construction and improvement of urban green space environment.

Nitrogen deposition caused higher increases in plant-derived organic carbon than microbial-derived organic carbon in forest soils.

Plant- and microbial-derived organic carbon, two components of the soil organic carbon (SOC) pool in terrestrial ecosystems, are regulated by increased atmospheric nitrogen (N) deposition. However, the spatial patterns and driving factors of the responses of plant- and microbial-derived SOC to N deposition in forests are not clear, which hinders our understanding of SOC sequestration. In this study, we explored the spatial patterns of plant- and microbial-derived SOC, and their responses to N addition and elucidated their underlying mechanisms in forest soils receiving N addition at four sites with various soil and climate conditions. Plant- and microbial-derived SOC were quantified using lignin phenols and amino sugars, respectively. N addition increased the total microbial residues by 20.5% on average ranging from 9.4% to 34.0% in temperate forests but not in tropical forests, and the increase was mainly derived from fungal residues. Lignin phenols increased more in temperate forests (average of 63.8%) than in tropical forests (average of 15.7%) following N addition. The ratio of total amino sugars to lignin phenols was higher in temperate forests than in tropical forests and decreased with N addition in temperate forests. N addition mainly regulated soil microbial residues by affecting pH, SOC, exchangeable Ca2+, gram-negative bacteria biomass, and the C:N ratio, while it mainly had indirect effects on lignin phenols by altering SOC, soil C:N ratio, and gram-negative bacteria biomass. Overall, our findings suggested that N deposition caused a greater increase in plant-derived SOC than in microbial-derived SOC and that plant-derived SOC would have a more important role in sequestering SOC under increasing N deposition in forest ecosystems, particularly in temperate forests.

Isoetes candelariensis, a new species of Isoetaceae (Lycopodiidae) from Argentina.

Isoetes candelariensis is a new species of Isoetaceae from Misiones, Argentina. This species is ephemeral, and grows on basaltic bedrock outcrop pools, in Urutau Reserve from Candelaria Department. This taxon differs from other aquatic Isoetes in this region by the unique combination of characters of its leaves, ligule, labium and megaspore ornamentation.

A new R package to parse plant species occurrence records into unique collection events efficiently reduces data redundancy.

Biodiversity data aggregators, such as Global Biodiversity Information Facility (GBIF) suffer from inflation of the number of occurrence records when data from different databases are merged but not fully reconciled. The ParseGBIF workflow is designed to parse duplicate GBIF species occurrence records into unique collection events (gatherings) and to optimise the quality of the spatial data associated with them. ParseGBIF provides tools to verify and standardize species scientific names according to the World Checklist of Vascular Plants taxonomic backbone, and to parse duplicate records into unique 'collection events', in the process compiling the most informative spatial data, where more than one duplicate is available, and providing crude estimates of taxonomic and spatial data quality. When GBIF occurrence records for a medium-sized vascular plant family, the Myrtaceae, were processed by ParseGBIF, the average number of records useful for spatial analysis increased by 180%. ParseGBIF could therefore be valuable in the evaluation of species' occurrences at the national scale in support for national biodiversity plans, identification of plant areas important for biodiversity, sample bias estimation to inform future sampling efforts, and to forecast species range shifts in response to global climate change.

A case study on the early stage of Pinus nigra invasion and its impact on species composition and pattern in Pannonic sand grassland.

Alien woody species are successful invaders, frequently used for afforestation in regions like semi-arid lands. Shrubs and trees create important microhabitats in arid areas. Understorey vegetation in these habitats has unique species composition and coexistence. However, the impact of solitary woody species on understorey vegetation is less understood. This study evaluated the effect of native (Juniperus communis) and invasive solitary conifers (Pinus nigra) on surrounding vegetation, where individuals were relatively isolated (referred to as solitary conifers). The field study conducted in Pannonic dry sand grassland in 2018 recorded plant and lichen species presence around six selected solitary conifers. Composition and pattern of understorey vegetation were assessed using 26 m belt transects with 520 units of 5 cm × 5 cm contiguous microquadrats. Compositional diversity (CD) and the number of realized species combinations (NRC) were calculated from the circular transects. Results showed native conifer J. communis created more complex, organized microhabitats compared to alien P. nigra. CD and NRC values were significantly higher under native conifers than invasive ones (p = 0.045 and p = 0.026, respectively). Native species also had more species with a homogeneous pattern than the alien species. Alien conifers negatively affected understorey vegetation composition and pattern: some species exhibited significant gaps and clusters of occurrences along the transects under P. nigra. Based on our study, the removal of invasive woody species is necessary to sustain habitat diversity.

Coordination of pinna, petiole, and root anatomical traits in 24 tropical-subtropical fern species.

Ferns are primitive vascular plants with diverse morphologies and structures. Plant anatomical traits and their linkages can reflect adaptation to the environment; however, these remain are still poorly understood in ferns. The main objective of this study was to explore whether there was structural coordination among and within organs in fern species. We measured 16 hydraulically related anatomical traits of pinnae, petioles, and roots of 24 representative fern species from the tropical and subtropical forest understory and analyzed trait correlation networks. In addition, we examined phylogenetic signals for the anatomical traits and analyzed co-evolutionary relationships. These results indicated that stomatal density and all petiole anatomical traits exhibited significant phylogenetic signals. Evolutionary correlations were observed between the tracheid diameter and wall thickness of the petiole and between the water transport capacity of the petiole and stomatal density. Conversely, anatomical traits of roots (e.g., root diameter) showed no phylogenetic signals and were not significantly correlated with those of the pinnae and petioles, indicating a lack of structural coordination between the below- and above-ground organs. Unlike angiosperms, vein density is unrelated to stomatal density or pinna thickness in ferns. As root diameter decreased, the cortex-to-stele diameter ratio decreased significantly (enhanced water absorption) in angiosperms but remained unchanged in ferns. These differences lead to different responses of ferns to climate change and improve our knowledge of the water adaptation strategies of ferns.

High preseason temperature variability drives convergence of xylem phenology in the Northern Hemisphere conifers.

Wood growth is key to understanding the feedback of forest ecosystems to the ongoing climate warming. An increase in spatial synchrony (i.e., coincident changes in distant populations) of spring phenology is one of the most prominent climate responses of forest trees. However, whether temperature variability contributes to an increase in the spatial synchrony of spring phenology and its underlying mechanisms remains largely unknown. Here, we analyzed an extensive dataset of xylem phenology observations of 20 conifer species from 75 sites over the Northern Hemisphere. Along the gradient of increase in temperature variability in the 75 sites, we observed a convergence in the onset of cell enlargement roughly toward the 5th of June, with a convergence in the onset of cell wall thickening toward the summer solstice. The increase in rainfall since the 5th of June is favorable for cell division and expansion, and as the most hours of sunlight are received around the summer solstice, it allows the optimization of carbon assimilation for cell wall thickening. Hence, the convergences can be considered as the result of matching xylem phenological activities to favorable conditions in regions with high temperature variability. Yet, forest trees relying on such consistent seasonal cues for xylem growth could constrain their ability to respond to climate warming, with consequences for the potential growing season length and, ultimately, forest productivity and survival in the future.

Adherence to plant-based dietary patterns and anthropometric indices among primary school girls in Kerman: A cross-sectional study.

The objective of the study was to evaluate the association between adherence to plant-based dietary patterns and anthropometric indices among primary school girls in Kerman. This cross-sectional study included 330 girls aged 6-12. A reliable and validated dish-based food frequency questionnaire was used to collect dietary data. Weight, height, and mid-arm circumference were measured, and z-score charts from the World Health Organization for girls between the ages of 5 and 19 were utilized. We used Satija et al. method to calculate plant-based diet index scores. After adjusting for age and energy, participants in the higher tertile of the overall plant-based diet index (PDI) had a lower height-for-age z score (HAZ) (0.38±0.12 vs. 0.59±0.12 kg; P = 0.033). Higher unhealthful plant-based diet index (uPDI) scores were significantly associated with decreased HAZ in models 2 (p = 0.028) and 3 (p = 0.035). A higher PDI score was associated with lower odds of being underweight and overweight, respectively (Model 1: P trend = 0.007, <0.001; Model 2: P trend = 0.010, 0.001). A significant inverse association was found between healthful plant-based diet index (hPDI) scores and underweight risk in the crude and adjusted models. (Model 1: P trend = 0.021; model 2: P trend = 0.018; Model 3: P trend = 0.031). Higher uPDI scores were associated with increased odds of being overweight and obese in all three models (Model 1: p trend<0.001; Model 2: p trend<0.001; Model 3: p trend = 0.001). We concluded that children who followed a plant-based diet had lower odds of being overweight and obese. Higher scores on the hPDI were linked to a decreased risk of being underweight, while higher scores on the uPDI were associated with an increased risk of being overweight and obese. The study suggests healthy plant-based diet may benefit children's weight and growth.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers.

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Contrasting physiological strategies explain heterogeneous responses to severe drought conditions within local populations of a widespread conifer.

Understanding how trees prioritize carbon gain at the cost of drought vulnerability under severe drought conditions is crucial for predicting which genetic groups and individuals will be resilient to future climate conditions. In this study, we investigated variations in growth, tree-ring anatomy as well as carbon and oxygen isotope ratios to assess the sensitivity and the xylem formation process in response to an episode of severe drought in 29 mature white spruce (Picea glauca [Moench] Voss) families grown in a common garden trial. During the drought episode, the majority of families displayed decreased growth and exhibited either sustained or increased intrinsic water-use efficiency (iWUE), which was largely influenced by reduced stomatal conductance as revealed by the dual carbon­oxygen isotope approach. Different water-use strategies were detected within white spruce populations in response to drought conditions. Our results revealed intraspecific variation in the prevailing physiological mechanisms underlying drought response within and among populations of Picea glauca. The presence of different genetic groups reflecting diverse water-use strategies within this largely-distributed conifer is likely to lessen the negative effects of drought and decrease the overall forest ecosystems' sensitivity to it.

Application of coconut shell activated carbon filter in vertical subsurface flow constructed wetland for enhanced multi-metal bioremediation and antioxidant response of Salvinia cucullate.

Coconut shell activated carbon (CNSAC) was applied as a filter layer in hybrid vertical subsurface flow constructed wetland (H-VSSF-CW), in order to enhance the multi-metal removal efficiency of the constructed wetland (CW) and to reduce heavy metal accumulation on Salvinia cucullata. Treatment P + AC, (having CNSAC filter layer), showed 32, 21 and 34% more Cd, Cr, and Pb removal efficiency than treatment P (without CNSAC layer). CNSAC activated carbon adsorbed Cd and Pb and Cr by functional groups -NH, -NO2, -C-O, -OH and -CO, and significantly reduced Cd and Pb exposure to S. cucullate. Chromium adsorption by CNSAC filter layer was half (just 25% of total input) of the Cd and Pb. In treatment P, due to high Cd, Pb and Cr accumulation in S. cucullate, the antioxidant defense mechanism of the plant was collapsed and cell death was observed, which in turn has resulted reduced biomass gain (5% reduction). On the other hand, in treatment P + AC, an antioxidant defense mechanism was active in the form significantly (p ≤ 0.05) increased of SOD, CAT and proline content while reduced MDA, EL, %EB and soluble sugar. So, the application of CNSAC increased the heavy metal removal efficiency of H-VSSF-CW by adsorption of a considerable share of heavy metal and hence, reduced the heavy metal load/exposure to S. cucullate.

Chionobium takahashii, gen. et sp. nov., associated with snow blight of conifers in Japan.

Gremmenia abietis (Dearn.) Crous (syn: Phacidium abietis) was originally described in North America to accommodate the species associated with snow blight of Abies and Pseudotsuga spp. In Japan, this species was first observed on the dead needles on Abies sachalinensis and Picea jezoensis var. jezoensis in 1969. However, the identity of Japanese species was unclear due to the lack of molecular data and the absence of anamorph description. In this study, we collected fresh specimens from various conifer species (A. sachalinensis, A. veitchii, Pic. jezoensis var. jezoensis, Pic. jezoensis var. hondoensis, Pinus koraiensis, and Pin. pumila) in Japan and revised the taxonomy based on morphological and phylogenetic analyses. Phylogenetic analyses based on nuc rDNA internal transcribed spacer ITS1-5.8S-ITS2 (ITS), nuc 28S rDNA (28S), and RNA polymerase II second largest subunit (RPB2) regions indicated that the species belongs to Phacidiaceae. Conidiomata formed in vitro produced pyriform, hyaline conidia without mucoid appendage, which distinguished the species from phylogenetically related genera. Consequently, we established Chionobium takahashii to accommodate the snow blight fungus in Japan. Further phylogenetic analyses also indicated that C. takahashii includes several distinct clades corresponding to the host genera (Abies, Picea, Pinus). Morphological differences among those clades were unclear, suggesting that C. takahashii may contain host-specific cryptic species.

Accumulation characteristics of liquid crystal monomers in plants: A multidimensional analysis.

Liquid crystal monomers (LCMs), identified as emerging contaminations, have been detected in soils and plants, but their accumulation characteristics in plants haven't been studied. Therefore, this study systematically investigated the accumulation characteristics of LCMs in plants from four dimensions (i.e., plant fruit species, soil types, plant growth stages, and LCMs categories) for the first time. The LCMs concentrations (9.96 × 10-4 to 114.608 ng/g) in 22 plant fruits were predicted by the partition-limited model. Grains with the highest lipid content showed the highest LCMs accumulation propensity. Plants grown in paddy soil showed a strong LCMs accumulation capacity. Results showed that the LCMs accumulation capacity in plants from soils decreased when the soil organic matter content increased. A preferential accumulation of LCMs in plant root systems during growth was found by the molecular dynamics simulations. Compared to polychlorinated biphenyls (as the reference contaminants of LCMs), LCMs exhibit higher accumulation in plant roots and lower translocation to shoots. For the fourth dimension, lipophilicity was found to be the main reason of LCMs accumulation by intergraded stepwise linear regression with sensitivity analysis. This is the inaugural research concentrating on LCMs accumulation in plants, providing insights and theoretical guidance for future LCMs management strategies multidimensionally.